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PREFACE 


The  present  investigations  form  a  part  of  a  study 
dealing  with  the  response  of  a  thin  spherical  shell  sub- 
ject to  an  internal  explosion,  being  carried  on  in  the 
Department  of  Engineering  Science  and  Mechanics  at  the 
University  of  Florida  under  a  research  contract  from  the 
U.S.  Army  Biological  Laboratories,  Fort  Detrick,  Maryland. 
The  general  problem  is  exceedingly  complex  and  requires 
effort  in  such  diverse  fields  as  detonation  chemistry, 
blast  waves,  stress  waves  in.  solid  media,  fracture,  vi- 
bration theory  and  reflection  phenomena.     It  has  been 
possible  to  deal  here  only  with  a  selected  number  of 
topics . 

Chapter  I  and  II  deal  with  the  calculation  of 
pressure  behind  the  shock  front  in  a  blast  wave  generated 
by  a  finite  source.     Theoretical  calculations  are  supple- 
mented by  experimental  data  obtained  at  the  University  of 
Florida.     In  Chapter  III,  the  nature  of  fracture  under 
impulsive  loading  is  examined.    An  empirical  criterion 
for  number  of  fragments  based  upon  the  experimental  data 


ii 


iii 


obtained  at  Fort  Detrick  is  proposed.     Finally,  in 
Chapter  IV,  pressure  calculations  are  related  to  the 
fracture  study  through  the  introduction  of  the  concept 
of  fracture  times. 

It  is  always  hard  for  a  graduate  student  to  list 
all  of  the  sources  of  his  intellectual  inspiration. 
Here,  only  some  more  pressing  obligations  can  be  met. 
During  two  years  of  work  on  this  problem,  the  author 
received  much  valuable  help  in  exploring  new  ideas  from 
Dr.  J.  Siekmann,   chairman  of  the  supervisory  committee. 
Mr.  D.  W.  Falconer  of  Army  Biological  Laboratories  was 
responsible  for  making  available  experimental  data  on 
fracture.     Mr.  D.  Mitchell  was  invaluable  in  conducting 
the  experimental  work  here.     Thanks  are  also  due  to 
Dr.  W.  A.  Nash,  Chairman  of  the  Department  of  Engineering 
Science  and  Mechanics  for  providing  financial  assistance. 
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NOMENCLATURE 


a  Radius  of  explosive  charge 

A  Dimensional  Constant 

B{&)  Constant  for  specified  Q 

Velocity  of  sound  in  air 

E  Energy  of  explosion 

f,  g,  h,  k           Nondimensional  flow  variables 

F,  G,  H,  K           Dimensional  flow  variables 

I  Impulse 

-l^  ,  m,  r\  '            Nondimensional  exponents 

M  Mass 

N  Number  of  fragments 

p  Pressure 

p  Dynamic  pressure 

r, 0  Spherical  Polar  coordinates 

R  Shock- radius 

t  Time 

v^,  u  Velocity  in  radial  direction 

,  v  Velocity  in  meridional  direction 

U  Velocity  of  shock  head 
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NOMENCLATURE  (Continued) 


Particle  velocity 

Critical  velocity  of  straining 

Ratio  of  specific  heats 

Similarity  variable 

=  (N/V) 

Angle  between  shock  radius  and 
normal  to  shock  head 


Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
in  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Doctor  of  Philosophy 
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Proceeding  from  general  mathematical  and  physical 
considerations,  a  group  of  similarity  transformations  has 
been  obtained  for  the  case  of  axisymmetric  flow  when  the 
Euler  equations  are  given  in  spherical  coordinates. 
Previous  similarity  solutions  due  to  G.  I.  Taylor  are 
shown  to  be  special  cases  of  these  transformations.  It 
is  also  shown  that  while  these  transformations  cannot 
be  applied  consistently  to  obtain  an  exact  solution  for 
the  case  of  a  blast  wave  generated  by  a  finite  source, 
they  may  be  used  to  find  an  approximate  solution. 
Numerical  results  are  presented  for  the  spherical  blast 
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wave  resulting  from  a  finite  source?  and  an  axisymmetric 
blast  wave  generated  by  a  point-source-finite-energy 
explosion. 

A  method  employing  a  mult i- spark  camera  device  for 
obtaining  the  Schlieren  record  is  presented.     Two  examples 
for  very  short  time  following  the  explosion  are  given. 
For  these  cases,  the  calculation  of  flow  variables  in  the 
contact  front  is  discussed. 

Rinehart  and  Pearson's  criterion  for  rupture  of 
cylindrical  shells  is  extended  to  spherical  shells.  An 
experimental  program  on  polystyrene  plastic  shells  of 
different  sizes  is  reported.     The  dominant  fracture 
pattern  observed  was  brittle.     Also,  critical  velocity  of 
straining  (V^j.)  was  found  to  be  directly  proportional  to 
the  cube  of  the  diameter  of  the  shell.     On  the  basis  of 

dimensional  considerations,  it  is  shown  that  shell  wall 
thickness  has  no  effect  on  the  number  of  fragments. 
Finally,  the  effect  of  material  properties  on  the  critical 
velocity  of  straining  and  the  number  of  fragments  is 
investigated. 

Fracture  time  is  defined  as  the  period  that  must 
elapse  between  the  impingement  of  the  blast  wave  on  the 
surface  of  the  shell  and  the  rupture  of  the  shell.  A 
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formula  for  calculation  of  fracture  time  is  obtained 
by  equating  the  impulse  imparted  to  the  final  momentum 
of  the  shell.     A  simplified  form  of  the  formula  is 
applied  for  calculation  of  fracture  times  for  thin  shells. 
It  is  observed  that  fracture  time  is  dependent  upon 
dynamic  pressure. 


CHAPTER  I 


APPROXIMATE  SOLUTION  OF  THE  PROPAGATION  OF  AN 
AXISYMMETRIC  BLAST  WAVE  GENERATED 
BY  A  FINITE  SPHERICAL  CHARGE 

1.  Introduction 

An  explosion  generally  consists  of  a  sudden 
release  of  a  finite  amount  of  energy  from  a  small  quanti- 
ty of  solid,   liquid  or  gas;  and  an  explosive  is  any  such 
solid,   liquid  or  gas.     Detonation  refers  to  the  propa- 
gation of  shock  wave  in  combustible  media,  whereas  blast 
refers  to  its  propagation  in  compressible  noncombustible 
media. 

Theoretical  and  experimental  investigations  have 
been  directed  on  the  one  hand  to  the  calculation  of  the 
equation  of  state  of  the  explosion  products  j^l,  2 J ,  and 
on  the  other  hand  to  the  propagation  of  shock  waves  in 
various  media.     Sir  G.  I.  Taylor     3 J  discussed  the  dynamics 
of  detonation  fronts  in  solid  explosives.     Recently  con- 
siderable interest  has  been  generated  in  detonation  in 
gases  due  to  its  applications  in  rocketry.     A  survey  of 
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progress  in  this  direction  is  contained  in  A.  Oppenheim, 
et  al.  IjiJ.     In  this  paper,  while  no  study  of  detonation 
as  such  is  undertaken,    it  is  demonstrated,   that  Taylor's 
solution  [3]  forms  a  special  case  of  the  transformation 
set  derived. 

( 

A  detailed  review  of  the  literature  concerning 
blast  waves  in  gases  was  given  by  I .   I.  Glass  [^5  J.  In 
the  following,   only  the  work  pertinent  to  the  present 
investigation  will  be  discussed.     G.  I.  Taylor  [^ej 
solved  the  problem  of  a  spherical  blast  wave  arising 
from  an  intense  point  source  explosion  in  air  and  pre- 
dicted that  pressure  will  be  attenuated  with  distance 
as  1/R-^,   where  R  denotes  the  shock  wave  radius.  While 
his  results  generally  agreed  with  observed  experimental 
phenomena  for  moderate  values  of  R,   agreement  was  not  as 
good  for  small  and  large  values  of  R.     The  same  problem 
has  also  been  solved  analytically  by  L.   I.  Sedov  [I'J, 
J.  von  Neumann  [sj ,   J.  L.  Taylor  [9]  and  R.  Latter  [loj. 
J.  A.  McFadden  [llj  tried  to  provide  solutions  for  short 
times  after  the  detonation.     His  method  is  applicable, 
however,   to  a  distance  of  only  five  per  cent  of  the 
explosive  radius.     H.  L.  Erode  |^12j  has  attempted  to 
provide  numerical  solutions  of  the  complete  equations  for 
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spherical  blast  waves.     Erode  predicts  that  pressure  will 
attenuate  at  a  rate  proportional  to  1/R  for  small  values 
of  R,   and  as  1/R^  for  moderate  values  of  R,  when  the 
explosion  products  have  absorbed  ambient  gases  approxi- 
mately ten  times  their  mass.     In  the  intermediate  range, 
pressure  decays  as  where  0.5   <^n   <Cl.5.  These 

results  for  a  particular  charge  of  T.N.T.  agree  reason- 
ably well  with  those  of  F.  J.  Berry  et  al.  fisj. 


formations  are  derived  for  axisymmetric  fluid  motion 
employing  Euler's  equations  in  spherical  coordinates, 
and  reducing  them  together  with  the  continuity  equation 
to  a  system  of  four  ordinary  differential  equations.  It 
is  shown  that  G.   I.  Taylor's  solutions  [^3,   6,   I4J  are 
special  cases  of  these.     The  applicability  of  the  Rankine- 
Hugoniot  conditions  at  the  shock  front  has  been  examined. 
It  is  demonstrated  that  it  is  possible  to  satisfy  them 
only  under  the  strong  shock  assumptions.  Numerical 
examples  are  provided  for  the  cases  of   (i)  spherical 
blast  waves  originating  from  a  finite  charge  and 
(ii)   axisymmetrical  blast  wave  generated  by  a  finite- 
energy-point-source  explosion. 


In  the  present  work,  general  similarity  trans- 
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2.     Equations  of  Motion  and  Boundary  Conditions 

The  blast  wave  is  assumed  to  originate  from  a 
sphere  filled  with  explosive  and  encircled  by  a 
diaphragm  around  its  periphery  in  the  equatorial  plane 
(see  Figure  1) .     It  is  suggested  by  the  geometry  that 
one  must  consider  the  variation  of  flow  variables  not 
only  in  the  radial  direction,  but  also  in  the    Q  ,  or 
meridional  direction.     Assuming  that  the  products  of 
explosion  can  be  approximated  by  a  perfect,  adiabatic, 
inviscid  gas,  the  equations  of  continuity,  momentum 
and  energy  may  be  expressed  in  spherical  coordinates 
as  follows: 


Equation  of  Continuity 


P  Cot 


Equations  of  Momentum 


4-  ^ 


4- 


iJ^  _  I// 


A    c^e  A 
1 


f  ^A 


(2) 


6 


5^  <Ll^  _j_  ^ 


Equation  of  Reversible  Adiabatic  Change  of  State 


(3) 


f 


r 


(4) 


l6  ^  ?^ 

where    lA  denotes  the  velocity  In  the  radial  direction, 
^&      the  velocity  in  the  meridional  direction,    \^  the 
pressure,    ^   the  density  and  V  the  ratio  of  specific 
heats.     In  real  gases,  the  specific  heat  of  a  gas  may  be 
approximated  by  a  power  series  in  temperature.     It  then 
becomes  possible  to  calculate  a  mean  specific  heat  for 
a  particular  temperature  range  from  the  relation 
(see  Reference  15) . 


1 


'g-'i  J 


r'Z 


Cf,  dT 


Using  this  relation,  a  mean  value  for  over  this 
temperature  range  can  be  determined 
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Thus,   for  air  in  the  range  of  temperatures  and  pressures 
usually  encountered  in  blasts  from  conventional  explo- 
sives =  1.2  gives  a  fairly  good  approximation. 

The  equations   (1-4)  are  to  be  solved  with  regard 
to  a  suitable  set  of  initial  and  boundary  conditions 
which  will  be  specified  later. 

Since  an  analytical  treatment  of  the  set  of  non- 
linear equations   (1-4)  presents  mathematical  difficulties 
beyond  the  present  power  of  analysis  one  has  to  rely  on 
numerical  or  similarity  methods.     A  powerful  technique 
is  finding  a  group  of  transformations  reducing  the  set 
of  partial  differential  equations  to  a  system  of  ordi- 
nary differential  equations.     In  what  follows,   the  trans- 
formations will  be  obtained  for  shock  waves  generated 
by  an  instantaneous  release  of  a  finite  amount  of  energy 
from  a  finite  source. 

A  first  step  in  the  process  of  seeking  similarity 
transformations  is  the  selection  of  a  suitable  similarity 
variable,   which  should  preferably  be  dimensionless . 
Since  the  flow  field  has  two  characteristics  lengths, 
i.e.,   the  Shockwave  radius  R  (     ,   t)  ,   and  the  radius 
locating  any  point  of  interest  in  the  flow  field,   r,  it 
would  then  seem  reasonable  to  assume  the  following  form 
of  similarity  variable: 
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A 


R  C6,t) 


The  next  step  is  to  seek  from  physical  reasoning, 
forms  suitable  for  the  dependent  variables.     In  the 
most  general  case  the  following  transformations  may 
be  used: 


c<*'    /3*'    T*  and  ^*  are  functions  of  r,    0  and  t. 
The  subscript  o  refers  to  the  ambient  atmosphere. 
However,  this  might  be  too  general  to  be  of  much  help  in 
a  specific  problem.     A  particular  subset  of  it  which  is 
useful  in  the  present  problem  may  be  written  as: 


-n 


) 


(5) 


where  n^^,  n^,  -£  and  m  are  positive  or  negative  constants. 
As  will  be  shown  below,   their  choice  is  not  completely 
arbitrary. 
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Substituting  from  (5)  into  the  continuity  equa- 
tion (1),  one  obtains 

o  u 

[ri  F'K  +riFK'i-  ZKFj-[^  ^  + 


R 


-h  KO  R  ^^'^    cot &=  0 


where  a  prime  denotes  differentiation  with  respect  to 
^  .     The  equation  can  be  made  independent  of  R  by 
putting 


(6) 


with 


(7) 


For  the  time  being,   nothing  more  will  be  said  about  the 
nature  of  A  and  B  (S) .     The  problem  of  their  evaluation 
will  be  treated  later. 

Another  useful  relation  between  ^,  m  and  n  is 
obtained  by  substituting  from  equations  (5-7)  into 
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the  first  equation  of  momentum  (2).     This  yields: 

where         -    V     P  /  f    denotes  the  velocity  of  sound 
o  o      o  o 

in  the  ambient  atmosphere.  To  make  this  equation  inde- 
pendent of  R,  the  following  is  required 

m  =  { -h  zn 

(8) 

No  additional  relation  between  and  n  can  be  obtained 
by  substitution  in  (3)  and  (4). 

From  (5-8) ,  the  general  set  of  transformations 
reducing  the  partial  differential  equations  of  flow 
(1-4)  for  the  flow  field  between  the  shock  wave  and  the 
origin  to  ordinary  ones  may  be  written  as  follows: 

u.  ^  Ti"^  F  in) )  v  =  v.'"^  G  (n) . 


(9) 

where  -C  and  n  are  arbitrary  constants. 
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Set  (9)  allows  a  double-infinite  manifold  of 
transformations  for  selected  values  of  /  and  n.  This 
point  may  further  be  illustrated  by  reference  to  the 
classical  Taylor  problem  of  spherical  detonation  and 
blast  waves. 

(i)  Spherical  detonation  wave  emerging  from 
geometric  center  of  the  explosive;     Assuming  the  deto- 
nation wave  to  be  spherically  symmetric,  one  has 
B  =  G  =  0  and  R  =  R(t)  .     Further    ^R/^  t  =  dR/dt  =  A  = 
constant.     Then  taking  n  =  /   =0,  equations   (9)  reduce 
to: 

dR  ^  /]  .  (10) 

dt 

The  velocity  of  the  detonation  wave  is  equal  to  the  time- 
derivative  of  R  i.e.  U  =  dR/dt  =  A.     Integration  of  the 
foregoing  equation  yields  R  =  At,   and,  consequently, 
^  =  r/At.     Introducing  now  a  new  variable    ^  =  A      =  r/t, 
equations   (10)  may  be  rewritten  as: 

These  expressions  were  employed  by  G.  I.  Taylor  ^3,  I4J 
to  solve  the  problems  of  the  air  wave  surrounding  an 
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expanding  sphere,  and  the  motion  of  products  of  an 
explosion  behind  spherical  detonation  waves  in  solid 
explosives . 

(ii)  Spherical  blast  wave  resulting  from  point- 
source;     Since  the  blast  wave  is  spherically  symmetric, 
one  can  put  B  =  G  =  0  and  R  =  R(t) .     Assuming        =  0 
and  n  =  3/2,   equations  (9)   reduce  to 


This  set  of  transformations  was  first  used  by 
G.  I.  Taylor  /sjto  reduce  the  system  of  partial  differ- 
ential equations  of  fluid  flow  to  a  system  of  ordinary 
ones . 


extremely  interesting  features.     Any  arbitrary  selection 
of  't.  and  n  based  upon  the  physical  nature  of  the  prob- 
lem yields  a  new  set  of  transformations.     A  still  more 
remarkable  property  of  equations   (9)   is  the  fact  that 
u  and  V  are  attenuated  with  R  as  the  shock  head  velocity 
a  result  upheld  and  suggested  by  intuition. 


(11) 


Examination  of  equations  (9)  reveals  several 
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Also,  if  t  is  taken  to  be  zero,  the  pressure  attenuates 
with  R  as  the  square  of  the  velocity  of  u  or  v. 

Attention  will  now  be  directed  to  transforming 
the  partial  differential  equations  (1-4)  to  ordinary 
differential  equations  for  the  problem  at  hand.  If 
/   =  0  the  substitution  of  (9)   in  (1-4)  yields: 

-firiK'/K+ZF/ri  +  F'+FK'/K-^j^l-QBriK'/K 
_  eriG' i-Q  cot&]=  0  ^^^^ 

-  A[n  F-r^lf'J  +  ^f^'-^^  L-n  BF-BqF'J 

- [n  +  tlc^'/G]-  [ft  ^  BG/n ]+F  [(^'/(^  >-  V^] 

''z  n  Bli  -i-Bq  h'1 


(14) 


+  ^[h'/H  -  r  k'/k]  =  0 


(15) 
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It  is  useful  at  this  point  to  make  the  equations 
dimensionless .     This  eliminates  A  from  the  equations 
(12-15) .     Introduce  the  following  nondimensional 
variables: 


A  ^  (16) 


Substitution  of   (16)   in  equations   (12-15)  yields: 

.rik'/k  +Z  i/ri  -t  f  V  f +  V?  ■[-3Bri 
-[nf-i-ri  f']^  ff'^  |_^_  nBf-B^i'J 

=  y^olk-J  i'lnsn  -hBrih'] 


(19) 


Cr^fe/zi  -zn-rih'/hj.[i  +8g/72j 


(20) 
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The  solution  of  equations  (17-20)  in  conformity  with, 
a  suitable  set  of  initial  and  boundary  conditions  should 
give  the  flow  field  between  the  origin  and  the  shock 
wave.     A  numerical  solution  of  these  equations  will  be 
the  next  objective  in  the  investigation.     However,  be- 
fore doing  this,  it  is  necessary  to  pause  to  examine 
the  question  of  initial  and  boundary  conditions. 

A  preliminary  examination  of  Euler's  equations 
(1-4)  reveals  that  one  must  require  for  t  7/  0:  ]/"  /     =  0, 
since  otherwise  the  last  term  in  equation  (1)  would  be- 
come infinite.     Further  the  normal  component  of  velocity 
along  the  diaphragm,   i.e.  v,  must  vanish.     More  precise- 
ly, for  t  7/  0  this  condition  can  be  formulated 
I/' I        ,  =0,  where  6         0  as  the  thickness  of 

diaphragm  approaches  zero.     In  the  limiting  case,  there- 
fore, this  condition  may  be  replaced  by      ^  jQ-TljZj  ~  ®* 
Also  u  and  v  must  vanish  at  r  =  0.     Thus,   the  natural 
boundary  conditions  for  u  and  v  are: 

L  0  =  rr/z 

(21) 
(22) 
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Integration  of  equations   (17-20)  can  be  carried 
out  numerically  step  by  step  if  for  any  value  of  the 
variable  ""^"i  the  functions  f,  g,  h  and  k  are  known. 
It  is  not,  however,   feasible  to  start  from  the  origin 

as  the  substitution  of  (21-22)  in  (17-20)  reveals  that 

f         I  f  / 

f(0),  g(0),  h(0)  and  k(0)  become  infinite  there.  There- 
fore, the  integration  must  be  started  from  a  point  other 
than  the  origin,   e.g.  the  shock  head. 

Any  specification  of  the  boundary  conditions  at 
the  shock  head  requires  a  knowledge  of  the  shock  surface 
shape.     It  seems  reasonable  to  assume  that  the  shock 
surface  is  a  function  of  both  S  and  t  and  that  R  =  R(^,t) 
gives  any  point  on  the  shock  surface.     Also,  then 
B  ((9)  =  1/R   <^R/<](9    may  be  calculated.     Next  the  proce- 
dure for  obtaining  values  of  f,  g,  h  and  k  at  r  =  R(6'  ,t) 
i.e.,        ^  1  ,  will  be  outlined. 

Now  the  Rankine-Hugoniot  conditions  at  the  shock 
front  (  ^  =  i.  )  are: 

r-i         z  ^  -^V^  ^  ^ 
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where  U  is  the  shock-front  speed,   and  subscripts   1_  and 
0  denote  the  quantities  iininediately  behind  and  in  front 
of  the  shock  wave  respectively.     However,  these  conditions 
cannot  be  satisfied  consistently  with  similarity  assump- 
tions |^6^.     Nevertheless,  these  may  be  replaced  by  the 
conditions  for  strong  shock  waves  when   ^1  / p' q  is  very 
large.  Thus, 

U^jlj  ^  Z/(t-t)-) 

where  X  =    n  /  \^  ^  and  u(l)  +  "J  v(l). 

• — ^  — ^ 

e^  and  e^     designate  unit  base  vectors  along  the  direc- 
tions of  increasing  r  and  Q   respectively.     These  rela- 
tions are  also  valid  when  Y"  is  replaced  by  Y"  (see 
Reference  15) . 

Under  the  general  transformations  (9),  these 
boundary  conditions  reduce  to  the  following: 

(23) 
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with 


Hence 


Jl^  COS  ?[  ,  s/M  ;^ 


(24) 


where  A  is  the  angle  between  the  radius  vector  and  the 
normal  to  the  shock  surface  R  =  R(  (9»   t) .     Note  that 

■k         It         it  -k 

f  ,  g  ,  h  ,  k    are  nondimensional  quantities,  but  not 
necessarily  the  same  as  f ,  g,  h,   and  k.     In  fact,  they 
are  related  to  dimensional  quantities  in  the  following 


manner: 


-t-zn 


(25) 


The  relations  (23)  represent  the  boundary  conditions 
for  the  equations  of  flow  which  have  been  transformed  to 
ordinary  differential  equations  under  the  generalized 
similarity  transformations   (9) .     It  is  seen  from  equa- 
tions  (23)   that  unless  ^  =  Q,  boundary  conditions  vary 
with  R,   thus  contradicting  the  similarity  hypothesis. 
Hence,   to  satisfy  the  similarity  hypothesis   ^  must  be 


zero. 
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This  is,  of  course,  what  was  assumed  earlier. 
When   -6=0,   equations   (25)  become  identical  with  equa- 
tions (16).     Hence  f * ,  g*,  h*  and  k*  become  also 
identical  with  the  corresponding  unstarred  variables. 
Consequently  from  equations   (23),   the  boundary  con- 
ditions for  (17-20)  are: 

frr  '''' 


where 

^=     f(i}  +  es^W  ,  iLoi^oj 


Before  proceeding  further,  a  few  remarks  concern- 
ing the  calculation  of  "A"  seem  to  be  in  order. 

The  total  energy  E  of  the  disturbance  consists  of 
two  parts,  the  kinetic  energy 

and  the  heat  energy 
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Since  it  is  simpler  to  consider  a  spherically  symmetric 
case,  attention  shall  be  devoted  in  the  following  pri- 
marily   to  it.     Of  course,  it  goes  without  saying  that 
the  remarks  below  will  apply  equally  to  the  axisymmetri- 
cal  case.     For  a  spherically  symmetric  shock-wave, 
equation  (27)   reduces  to: 


(28) 

From  this  relation  it  is  clear  that  unless  71  =  3/2,  E 
depends  explicity  on  R{E  c?(  R3-2nj  hence  through  it 

implicitly  on  time  t.     The  case  Vi  =  3/2  corresponds  to 
the  point  source  finite  energy  explosion. 

The  problem  that  now  faces  one  is  "what  happens 
when  Tl  ^  3/2  and  lies  between  0  ^  ?^  T3/2  i.e.  for 
short  moments  immediately  after  explosion? " 

E  must  equal  the  energy  of  chemical  explosion,  and 
at  the  same  moment  must  depend  upon  t  if  similarity  rela- 
tions are  to  be  satisfied.     This  point  has  also  been 
discussed  by  H.  Morgenroth  \_}-^\  and  M.  H.  Rogers  |^17^. 
At  this  point,  it  is  instructive  to  examine  the  physical 
nature  of  the  blast  wave  generated  by  a  finite  source. 
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Flow  behind  the  leading  shock  in  this  case  is  not 
continuous  [l3.  is]. 

The  explosive  gases  behave  like  a  piston  and  feed 
energy  into  the  eimbient  gases  behind  the  main  shock- 
front.     The  piston  surface  is  a  contact  surface.  Be- 
sides the  contact  surface,  there  must  occur  additional 
shocks  as  shown  in  Figure  2.     Initially,  contact  surface 
and  leading  shock  are  together  and  region  I  is  non- 
existent.    However,   later  on  as  ambient  gases  are  absorb- 
ed by  the  blast  wave,   region  I  develops. 

It  is  essential  to  find,  therefore,  three  separate 
solutions.     The  method  outlined  here  can  only  be  used  to 
find  solution  in  region  I.     The  solution  cannot  be  extend- 
ed beyond  the  contact  surface  and  hence  it  is  not  possible 
to  satisfy  the  inner  boundary  condition  of  zero  velocity 
at  the  center.     Thus  in  region  I,   the  relation  E^o(^B?~^^ 
is  satisfied.     Theoretically,  then,  by  changing  limits  of 
integration  in  equation  (28)  so  as  to  include  only  region 
I,   it  should  be  possible  to  calculate  A.     This,  however, 
requires  finding  solutions  for  regions  II  and  III  as  well. 
The  energy  in  all  three  regions  equals  total  energy  (E) , 
Since,  however,   it  is  not  possible  to  apply  similarity 
solution  to  regions  II  and  III,   it  is  now  necessary  to 
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t 

Figure  2:     Flow  behind  the  Main  Shock 
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abandon  the  equation  of  energy  to  determine  "A".     A  must 
be  found  from  experimental  data  concerning  the  time 
history  of  shock  propagation. 

3.     Numerical  Calculations  for  Spherical  Blast  Waves 

In  this  section  it  is  proposed  to  provide  a  solu- 
tion for  the  spherical  blast  waves  resulting  from  both 
the  finite  source  and  point  source  explosions.  In 
Section  I  certain  previous  theoretical  as  well  as  experi- 
mental results  were  reviewed.     On  the  basis  of  these,  it 
is  possible  to  apply  different  similarity  solutions  for 
different  times. 

The  equations  of  motion  in  dimensionless  form  are 
obtained  from  equations   (17-20)  by  ignoring  the  Q 
coordinate.     These  may  be  written  after  some  elementary 
algebraical  manipulations  as: 


rlll:fl^_  iLi^  f  rI]J  ^n-^z- 
_  2  f -ZnrilY- 


(30) 


2nk__ 

r(n'f) 


(31) 
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where    Y'  and  y   denote  the  values  of  the  ratio  of 

o 

specific  heats  behind  and  in  front  of  the  shock  wave, 
respectively. 

It  will  be  assumed  in  the  following  that  for 
small  times  the  pressure  behind  the  shock  front  is  con- 
siderably greater  than  in  front  of  it  in  order  to 
justify  the  use  of  strong  shock  boundary  conditions. 
From  equation  (26),   the  boundary  conditions  at  ^  =  1 
are: 

(32) 

The  solution  of  equations   (29-31)   subject  to  the  con- 
ditions  (32)  can  be  carried  out  only  by  numerical  means. 

The  numerical  method  employed  to  solve  equations 
(29-31)   is  a  special  form  of  the  Runge-Kutta  method,  due 
to  Merson  [^^^J-     The  differential  equations  (29-31)  were 
programmed  on  the  IBM  709  computer,  utilizing  the  numeri- 
cal method  mentioned  above  to  yield  results  for  various 
values  of  n,  within  three  per  cent  accuracy  at  each  step. 
The  integration  interval  employed  was  w  =  0.002  (except 
in  Case  3) .     Results  have  been  obtained  for  the  following 
cases: 
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Case  1 


n  -  0-5 


Behavior  of  compressed  ambient  gas  region  for  short 
times  after  explosion . 

Case  2 


Same  as  in  Case  1  expect  in  so  far  as  it  relates  to 
intermediate  values  of  time. 

Case  3 


The  integration  interval  employed  was  0.01.     This  is 
identical  with  Taylor's  finite-energy  point  source 
explosion  solution.     In  Table  1,  the  solution  given  by 
Taylor      J  is  compared  to  the  one  obtained  here.  The 
other  results  of  these  computations  are  reproduced  in 
Figures  3-5. 

If  one  regards  the  contact  surface  as  the  outer 
boundary  of  an  expanding  sphere  of  radius   (r) ,  then 
its  velocity  is  given  by  (see  Reference  14) : 


T^To  ■=  l-"^  )  1-0 


n:=  1-5 


26 


TABLE  1 


SPHERICAL  SYMMETRIC  POINT  SOURCE  EXPLOSION: 
COMPARISON  OF  TAYLOR'S   [s]  RESULTS 
WITH  THE  AUTHOR'S 


h 

f 

k 

V 

i 

Taylor 
Solution 

Present 
Solution 

Taylor 

Present 

Taylor 

Present 

1.00 

1.167 

1.167 

0.833 

0.833 

6.000 

6.000 

0.98 

0.949 

0.955 

0.798 

0.798 

4.000 

4.068 

0.96 

0.808 

0.813 

0.767 

0.766 

2.808 

2.887 

0.94 

0.711 

0.714 

0.737 

0.736 

2.052 

2.122 

0.92 

0.643 

0.644 

0.711 

0.709 

1.  534 

1.603 

0.90 

0.593 

0.  592 

0.687 

0.684 

1.177 

1.235 

0.88 

0.556 

0.553 

0.665 

0.661 

0.919 

0.967 

0.86 

0.528 

0.524 

0.644 

0.639 

0.727 

0.765 

0.84 

0.507 

0.  502 

0.625 

0.611 

0.578 

0.619 

0.82 

0.491 

0.484 

0.607 

0.600 

0.462 

0.490 

0.80 

0.478 

0.471 

0.590 

0.582 

0.370 

0.  395 

0.78 

0.468 

0.460 

0.573 

0.565 

0.297 

0.319 

0.76 

0.461 

0.452 

0.  557 

0.  548 

0.239 

0.257 

0.74 

0.455 

0.446 

0.  542 

0.532 

0.191 

0.208 

0.72 

0.450 

0.441 

0.527 

0.516 

0.152 

0.168 

0.70 

0.447 

0.437 

0.513 

0.501 

0.120 

0.135 

0.68 

0.444 

0.434 

0.498 

0.485 

0.095 

0.108 

0.66 

0.442 

0.432 

0.484 

0.463 

0.074 

0.086 

0.64 

0.440 

0.430 

0.470 

0.455 

0.058 

0.068 

0.62 

0.439 

0.429 

0.456 

0.441 

0.044 

0.054 

0.60 

0.438 

0.428 

0.443 

0.426 

0.034 

0.042 

0.58 

0.438 

0.427 

0.428 

0.419 

0.026 

0.033 

0.56 

0.437 

0.426 

0.415   !  0.396 

0.019 

0.025 

0.  54 

0.437 

0.426 

0.402 

0.382 

0.014 

0.019 

0.52 

0.437 

0.425 

0.389 

0.367 

0.010 

0.015 

0.50 

0.436 

0.425 

0.375 

0.352 

0.007 

0.011 

0.45 

0.425 

0.315 

0.005 

0.40 

0.425 

0.278 

0.002 

0.35 

0.425 

0.240 

0.001 

Q.30 

0.425 

0.201 

0.000 

0.25 

1  1 

0.425 

0.159 

0.000 

0.20 

0.425 

0.132 

0.000 

0.15 

j  1 

0.425 

0.053 



0.000 

0.12 

 f 

0.425 

0.001 

0.000 

0.00* 

 1 

0.425 

0.000 

0.000 

This  entry  was  extrapolated. 
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Figure  5:     Plot  of  k  versus  1^ 
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Also  from  equations   (9)  and  (16) , 


u.^  u  {(n) 


Comparing  the  two  relationships,  the  contact  front  is 


Thus  for  the  Cases  1  and  2,  the  contact  front  is  speci- 
fied by 


Evidently,  any  solution  beyond  these  points  is  not 
physically  realistic. 

It  was  shown  previously  that  for  finite  source 
explosions,  the  calculation  of  "A"  must  be  made  from 
experimental  results.     Hence,  it  is  not  possible  here  to 
cite  any  results  for  Cases  1  and  2.     "A"  can  be  determined 
in  these  cases  in  the  following  manner.     From  Schlieren 
photographs,   it  is  possible  to  obtain  the  velocity  of 
propagation  U  of  the  main  shock.     Once  U  is  known,  A  may 
be  obtained  by  application  of  the  formula  U  =  AR"^. 
Experimental  investigations  along  these  lines  are  report- 
ed in  Chapter  II. 


defined  by 


Case  1 


Case  2 


3X 


The  situation  is  different  with  regard  to  Case  3. 
Here,   equation  (28)  reduces  to  the  following  simple  form: 


Numerical  integration  yields 

Compare  this  to  the  values  obtained  by  G.  I.  Taylor  [^6]. 

His  computations  show  that  E  =  5.36    T        and  thus  are 

^o 

in  quite  remarkable  agreement  with  present  results. 

D.  L.  Jones  [^2oJ  gives  a  slightly  different  value  (E  =  5.33 

n  2 

J        ) .     Finally,  given  E  and    f^,  A  may  be  evaluated. 

4.     Numerical  Calculations  for  Axisymmetrical  Point- 
Source-Finite-Enerqy  Blast  Wave 

The  equations  of  motion  (17-20)  may  be  written  in 

a  form  more  suitable  for  numerical  procedure  as  follows: 

-2f  +  8f^/«r  _  _a!_  2nh  -t-  gV^ 


(33) 


(34) 
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71  6- 

where 


(35) 


(36) 


Also,  from  equations  (24)  and  (26),  the  boundary  condi- 
tions at  ^  =  1  are: 

f  -  Cos  A     ■)  3  =  —  A 

h  -  ,     h  -    l-^^  (37) 

^  ~  i-^r  ^  r-i 

The  numerical  technique  employed  in  the  programming 
of  equations   (33-36)   for  solution  on  IBM  709  Computer 
was  the  same  as  described  in  Section  3.     Also,   the  same 
integration  interval  w  =  0.002  was  employed. 

It  is  assumed  that  the  explosion  is  due  to  a  point 

source  in  air.     Therefore,         ,    f    and  AZ  have  the 

'  o 

values  1.200,   1.200  and  1.500,  respectively. 
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To  preform  numerical  calculations,  it  is  necessary 
to  specify  R{B  ,  t)  at  some  moment  so  as  to  make  the 
calculation  of  B(^  )  and  ^  feasible.     It  may  happen  in 
an  actual  problem  that  B  depends  on  t  as  well.     In  this 
case,  it  will  be  necessary  to  specify  the  shape  of  the 
shock  wave  R(^,  t)  throughout  the  time  interval.  The 
similarity  solution  may  be  applied  now  by  approximating 
B  in  a  short  time  interval  /St,   so  that  B  depends 
essentially  on  3  only.     This,  though  interesting  will 
not  be  pursued  further  in  the  following. 

While  the  exact  shape  of  the  shock  head  must  be 
determined  experimentally,  it  was  thought  fit  to  exeimine 
here  the  hypothetical  case  when  the  shock  head  happens 
to  be  an  ellipsoid  of  revolution.     This  can  be  done  if 
one  assumes  that  the  ellipsoid  of  revolution  approxi- 
mates the  shock  surface  in  a  short-time  interval  A  t.  A 
cross- sectional  view  is  given  in  Figure  6.     The  equation 
of  an  ellipse  can  be  written  in  the  form 

(38) 

where  R  =  R{  B ,   t)  denotes  the  radius  of  the  shock  head, 
the  semi-major  axis,  and  Rq        the  semi-minor  axis. 
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Clearly  and  Rq  -  ^  are  functions  of  time.  Then,  from 
equations  (38) 

Q^d^yR  ^  s/A/  e  cos  e Ocno-S)^~_Rf  J 

\  =     /^RC  TW  (-8)  -^-B  ^''^ 

The  calculations  were  conducted  for  the  special 
case  when        =  2.2  and  %  -  %       2  (Units  are  not 
important.     It  is  only  the  ratio  Rq/Ro"^  ^^^^  matters.) 
Then  B(^),   f(l)  and  g(l)  can  be  calculated  from  equa- 
tions  (39).     Below  are  tabulated  some  of  these  values: 

TABLE  2 


B, 

f  AND  g  AS 

FUNCTIONS  OF 

& 

^o 

-  Y  =  1.1 

B 

B(6>) 

fd) 

0°,  90° 

-0.000 

0.9091 

0.0000 

30° 

-0.162 

0.9057 

0.0783 

60° 

-0.154 

0.9063 

0.0712 

The  differential  equations  were  solved  for  these  particu- 
lar values  of  S  .     Results  are  reproduced  in  Table  3. 
Some  explanatory  notes  are  required  to  fully  comprehend 
Table  3.     These  comments  may  be  divided  into  two  parts, 
i.e.,  concerning  calculations  for 
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0  =   O',  90' 


(a) 

and 

O     =  qQ.  90° 
In  this  case,  B  =  g  =  0  and  consequently  the 
system  of  equations   (33-36)  becomes  identical  with  that 
for  spherically  symmetrical  shock  waves.  Calculations 
can  be  carried  on  to  rather  low  values  of  ^  (;i:;0.1),  when 
the  accumulated  error  necessitates  approximation.  Since 
this  has  already  been  covered  in  Section  3.  nothing  more 
need  be  said  here. 

0   =  30°,  60Q 
The  numerical  calculations  in  this  case  are  to  be 
carried  out  in  three  essential  steps.     First,  the  system 
(33-36)  is  examined.     The  solution  can  be  carried  out  up 
to    ^  =  0.974,  when  the  small  value  of  g  makes  it 
necessary  to  change  the  numerical  procedure  (if  this  is 
not  done,   there  results  an  error  in  the  value  of  g  of 
its  own  order) .     m  the  second  step,  g  is  put  equal  to 
zero.     This  means  equation   (35)  drops  out  and  system 
(33-36)  becomes  identical  with  that  for  the  symmetrical 
case.     This  approximation  allows  the  solution  up  to  the 
point  0.6.     At  this  point,  however,  k  :^  0  and  h 
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approaches  a  constant  value,     f    must  approach  zero  at 

"7    =0.     Between     ''^  ~  0.6  and  ^  =  o,   f   >^  ^  provides 
a  good  approximation.    Also,  around   TJ^^    0.6,  k*  be- 
comes small  such  that  it  may  be  neglected.     Then  the 
approximation  f  =      ,  yields 

h'=  3/2  r^k 

The  calculation  of  A  is  straightforward  in  this  case 
as  E  =  const  (Section  2) .     Evaluation  of  the  integral 
(27)  leads  to  the  result 

E  =  11.^37  fo 

Hence  given  the  values  of    /    and  E,  A  may  be  found. 

o 

Calculations  have  also  been  carried  out  for  the  case 
when    7^=  = 
in  Table  4.  Here 

E=  ^-6  91  f^fl^ 

5 .  Conclusions 

The  existence  of  a  general  set  of  similarity  trans- 
formation for  the  axi symmetrical  Euler  equations  (in 
spherical  coordinates)  has  been  shown.     Previous  simi- 
larity solutions  of  Taylor  /I,  6,   I4J  have  been  demon- 
strated to  be  special  cases  of  these  transformations.  It 
is  not  possible  to  satisfy  the  exact  Rankine-Hugoniot  con- 
ditions at  the  shock-front.     However,  the  conditions 


Q      -.4  and       =  1.5.     These  are  reproduced 
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at  the  shock  front  can  be  satisfied  under  strong  shock 
assumptions . 

When  the  source  of  blast  wave  is  finite,   it  is 
not  possible  to  use  the  energy  relationship  to  solve  the 
problem  completely.     Despite  this  difficulty,   the  method 
has  been  demonstrated  to  be  capable  of  providing  approxi- 
mate solutions  in  such  cases,  when  information  on  time- 
history  of  the  shock  front  is  available.     Finally,  it 
is  possible  to  apply  this  method  to  analyze  the  problems 
of  axisymmetric  blast  waves.     In  this  instance,  however, 
one  must  require  additional  information  about  the  shape 
of  the  shock-head. 


CHAPTER  II 


SOME  EXPERIMENTAL  RESULTS  FOR  BLAST 
WAVES  ORIGINATING  FROM  A 
FINITE  SOURCE 

1.  Introduction 

In  Chapter  I  a  similarity  solution  has  been 
obtained  for  the  compressed  air  region  between  the 
main  shock  and  the  contact  surface  in  a  spherically 
symmetric  blast  wave  generated  by  a  finite  charge.  It 
was  pointed  out  that  in  order  to  complete  the  solution, 
it  is  necessary  to  determine  a  certain  constant  "A" 
from  experimental  data.     Once  "A"  is  known,   it  is  possible 
to  find  the  pressure,  density  and  velocity  distribution 
behind  the  main  shock.     The  most  convenient  way  to  evaluate 
"A"  is  to  first  find  "U"  -  the  main  shock  velocity  - 
from  Schlieren  pictures,  and  then  apply  the  relation 
U  =  AR       where  R  is  the  shock  wave  radius  and  n  is  a 
constant.     In  this  Chapter  a  method  employing  a  multi- 
spark  camera  for  obtaining  the  Schlieren  record  is  pre- 
sented.    Also,   two  examples  for  very  short-time  following 
the  explosion  are  given. 
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2.     Equipment  and  Procedure 

The  multi-source  spark  camera  is  of  the  so-called 
Cranz-Schardin  type  |j2lj,   and  was  first  developed  at  the 
National  Physical  Laboratory  (NPL) ,  England  [22] . 
Reference  [^23^  describes  a  modified  version  of  the  NPL 
model,  being  used  at  the  University  of  Toronto  Institute 
of  Aerospace  Studies,  Canada,   for  visualizing  flow  be- 
hind a  blast  wave.     The  camera  and  the  auxiliary  system 
(Figure  7)  employed  at  the  University  of  Florida  is 
basically  the  same  as  the  one  used  at  UTIAS. 

The  camera  essentially  consists  of  a  number  of 
individual  Schlieren  systems.     It  makes  use  of  a  pair  of 
common  main  mirrors  and  groups  of  distinct  but  identical 
light  sources,   lenses  and  knife  edges.     Each  separate 
Schlieren  system  is  triggered  at  a  predetermined  time 
from  an  initial  pulse  originating  from  the  flow  phenome- 
nom  to  be  observed.     In  the  present  case,  a  given  flow 
is  visualized  from  five  slightly  different  angles  and 
recorded  as  five  different  distinct  photographic  Images 
on  a  single  sheet  of  film.     The  time  between  exposures 
can  be  varied  by  means  of  short-duration  pulse  techniques 
and  time-delay  units,   from  a  minimum  of  2  ^seconds 
upwards.     The  chief  disadvantage  of  this  camera  is  the 


Figure  7:     Test  Stand 
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limited  number  of  pictures.     However,  as  pointed  out  in 
Reference    23  ,  this  disadvantage  is  offset  by  the 
simplicity  of  design,  economy,  resolution,  pre-set 
variable  or  constant  framing  rate  to  suit  a  given 
experiment. 

Two  different  sizes  of  uncased  PBRDX  charge  were 
employed  in  the  testing  program  i.e.,  (1)  a  spherical 
charge  of  3/8  inch  diameter  varying  in  weight  from  0.56 
grams  to  0.62  grams  and  (2)   a  spherical  charge  of  1/2 
inch  diameter  varying  in  weight  from  1.27  grams  to  1.31 
grams.     These  charges  were  detonated  by  use  of  a  high- 
energy  blast  cap  activated  by  a  one  microfarad  5000  volt 
condenser. 

The  time  delay  unit  was  set  to  trigger  the  first 
spark  source  after  4.3  microseconds  of  the  initiation  of 
the  blast  cap.     Succeeding  frames  were  taken  at  7,  10, 
15  and  20  microseconds,  respectively.     This  sequence  was 
adhered  to  throughout  the  entire  experimental  program  to 
obtain  reproducible  results.     Figure  8  reproduces  a 
typical  firing. 


Figure  8:     A  Typical  Firing  Sequence  for 
0.6  gm  PBRDX  Charge 
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3.     Experimental  Results 

A  close  look  at  Figure  8  shows  that  for  very  short 
times  immediately  following  the  explosion  such  as  in  the 
present  case,   the  second  shock  has  not  still  separated 
from  the  contact  surface.     There  is  a  very  thin  compressed 
air  region.     This  phenomenon  was  previously  observed  by 


H.  Schardin  |^24j .     Another  remarkable  feature  is  the 
irregularity  of  the  outer  boundary  of  the  expanding 
explosive  gases.     The  flat  region  from  which  two  wires 
protrude  is  due  to  the  effect  of  the  nonsymmetrical 
explosion  of  the  blast  cap. 

To  find  the  shock  wave  radius  R  as  a  function  of 
time  t,   it  was  found  necessary  to  approximate  the  irregu- 
lar material  boundary  by  a  circle.     Also,   the  flat  region 
was  neglected.     R  was  read  to  within  0.5  mm.    (scale  3.60  c: 
=  2  inches)  from  the  pictures.     The  values  for  R  are  given 
in  columns  2-4   (0.6  gm  charge)  and  columns  6-8  (1.3  gm 
charge)  of  Table  5.     Then  the  velocity  U  =  (dR/dt)  of  the 
shock-head  was  determined  by  approximating  the  differ- 
ential quotient  dR/dt  by  the  difference  quotient: 
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These  values  of  U  are  plotted  (see  Figure  9)  on  a 

log- log  paper.    The  experimental  points  show  an  even  scatter 
around  a  straight  line.     If  one  now  applies  the  formula; 
U  =  AR~^,   the  values  of  A  and  n  may  be  directly  read  from 
the  graph.     (Regarding  the  R-axis  as  the  abscissa  and 
the  U-axis  as  the  ordinates,   n  is  simply  the  slope  of  the 
straight  line.     The  value  of  A  equals  U  at  R  =  1  inch.) 

In  order  to  check  the  accuracy  of  the  straight- line 
approximation,   R  was  calculated  by  using  the  formula 

(this  formula  can  be  easily  obtained  by  integrating  the 
equation 

within  the  limits  t^,   R^  and  t,   R)  where  t^  and  R^ 
correspond  to  the  first  frame.     The  calculated  values 
of  R  are  given  in  columns  5   (0.6  gm  charge)  and  9   (1.3  gm 
charge)  of  Table  5.     They  differ  by  less  than  five  per 
cent  from  observed  values  of  R  -  an  extremely  good  agree- 
ment in  view  of  the  inaccuracies  in  reading  the  values  of 
R  from  the  pictures.     The  values  for  U  can  be  extrapolated 
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Figure  9:     Variation  of  Shock  Head  Velocity 
with  Shock  Radius 
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from  the  experimental  curves  by  simply  extending  the 
curves.     Care  must  be  taken  in  doing  this.     U  can  be 
approximated  in  this  manner  to  approximately  R  =  0.5 
inch  on  the  lower  side.    For  values  smaller  than 
R  =  0.50  inch,  the  approximation  is  bound  to  be  unreli- 
able as  at  this  point  the  shock-head  velocity  U  becomes 
equal  for  both  small  and  large  charges.     Similarly  on 
the  upper  end,   it  does  not  seem  possible  to  extrapolate 
beyond  R  =:  2.50  inches  due  to  low  overpressures. 

4 .     Expansion  of  Explosion  Products 

The  treatment  here  follows  closely  that  outlined 
by  Zeldovich  and  Kompaneets  [^2 J.     Given  the  initial  load- 
ing density         of  the  solid  explosive,  it  is  possible  to  ■ 
find  the  detonation  velocity  D  and  the  density  of  the 

explosion  products  in  the  detonation  front.     For  the  pre- 

jp  3 

sent  explosive,  with  a  loading  density  of  1.50  gm/cra  , 

the  detonation  velocity  D,  the  pressure  p^  in  the  detona- 
tion front,  and  the  density         have  the  values  8000  m/sec, 

3  3 
268.10    atm.  and  2,17  gm/cm    respectively.     The  explosion 

products  expand  in  air  in  essentially  two  steps: 

(i)  from  (  f  J  to  (        ,    f  )  according  to 

d        d     ^         '  X  X 

-  7^ 

polytropic  gas  law       ^  ^      1  =  const  and 
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(ii)  from  {   p  .    ?  )  to  any  lower  (   p.  f  ) 
according  to  poly tropic  equation  f       2  _  const. 

(  ^        •    ^il^         of  gives  appropriate 

expansion  in  earlier  stages,  its  use  after  a  certain 

point  tends  to  give  too  large  an  expansion.  This 

necessitates  the  use  of  ^1  (  ^  V^)   to  describe  the 

2       ^  1 

latter  stage.     The  juncture  point  represents  the 
boundary  line  between  the  two  stages.     The  juncture 
point  {  P  ,    ^  )  which  is  of  special  interest  in  the 

X  X 

expansion  of  the  explosion  products  in  air  if  found  by 
solving  the  equations 


and 


simultaneously.     The  heat  of  explosion  is  denoted  by  Q. 
In  the  present  case,   taking  Q  =  1280  cal/gm,    /"^^  =  3  and 

=  1«3,    (these  values  are  suggested  by  the  given 
quantities  and  D.     Also  see  Reference  2)  the 

Juncture  point  (  p  ,     f  )  is  found  to  be 

l^x  ^  1  5  S  0  czfm 
=   0-3  9  ^?v/c/n^ 
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5.    Pressure  Calculations 

From  Chapter  I,  pressure  p,  density  ^  and  velocity 
u  behind  the  main  shock  in  the  conpressed  air  region  are 
given  by: 


11=  uf(n^ 

where  p^,  ^  ^  refer  to  values  of  pressure  and  density  in 
the  ambient  atmosphere.     In  these  equations,  ^  =  r/R 
denotes  a  dimensionless  similarity  variable,  where  r 
locates  any  point  in  the  flow  field  from  the  center  of 
the  explosion. 

For  ^=1   (i.e.,   just  behind  the  main-shock)  under 
strong  shock  assumptions  h,  k  and  f  have  the  values: 

For  the  case,   r=  1.2  I.e.,  when  the  equation  of 
state  of  air  is  approximated  by  the  law  p  =  const, 

these  values  become: 

h(l)  =  I'O  91  ;  kU)=  11-000  ^  -f[l)^  0-009 


53 


It  is  possible  to  find  h(^),  k(^)  and  f{^)  in 
the  compressed  air  region  behind  the  main  shock  as  shown 
in  Chapter  I.     Also,  the  contact  front  may  be  determined. 
The  contact  front  is  thus  specified  by: 

(a)  0-/^         Chsr^e^  n=o.89 

/2  (^ O'  0  23 

(b)  1.3  C^^>^^     97-  0-70 

-f  C^)  ^  ^-  9^0^  o-  831 

Cq)  =  0-  349 

Across  a  contact  front,  pressure  (p)  and  velocity  (u)  are 
continuous  whereas  density  (^)  varies  discontinuously . 
The  density  across  the  contact  front  may,  then,  be  de- 
termined from  the  relation 

/-^  J  =  constant 

For  example  at  R  =  1.25  inches,  the  pressure  and  density 
across  the  contact  front  are  given  by 

(a)  0  ■  l£  Ch^r^ 

0-9  20:   fD=  (,7  a-/-^  f=0'(:)^^l^ 

(b)  1-3   f^ri.  cZ/i^y^e 

■  l:y^lO0Cfi7r2j=O'O^7f^ 

Similar  calculations  can  be  performed  for  other  values  of  R. 


CHAPTER  III 


ON  FRACTURE  OF  A  THIN  SPHERICAL  SHELL 
UNDER  BLAST  LOADING 

1.  Introduction 

The  behavior  of  materials  under  impulsive  loading 
is  of  great  interest  in  a  number  of  situations.  Apart 
from  the  military  applications   (e.g.  construction  of 
bombs) ,   it  is  important  in  the  understanding  of  various 
phenomena  involving  short-time  collisions   (e.g.  meteorite 
collisions) .     It  has  been  recognized  for  a  long  time  that 
materials  under  impulsive  loading  behave  in  an  essentially 
different  manner  than  in  the  static  situation.     In  the 
present  work,  an  experimental  investigation  was  directed 
towards  the  study  of  the  fracture  pattern  and  applica- 
bility of  a  modified  form  of  Rinehart  and  Pearson's 
criteria.     Experiments  on  thin  shells*  made  of  polystyrene 
plastic  revealed  the  fracture  pattern  to  be  primarily  of 

*Thin  shell  is  defined  as  the  one  for  which  wall 
thickness  is  less  than  10  per  cent  of  the  shell  radius. 
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a  brittle  type.     Also,  the  critical  velocity  of  strain- 
ing (V    )  was  found  to  be  directly  proportional  to  the 
cr 

cube  of  the  shell  diameter.     Finally  the  effect  of  shell 
wall  thickness  and  material  properties  on  the  number  of 
fragments  is  examined  on  the  basis  of  dimensional  analysis. 

2.     Fracture  Criteria 

The  first  mathematical  treatment  of  fracture  of 
elastic  media  was  given  by  A.  A.  Griffith  [^25jl,  and  was 
later  extended  independently  to  plastic  materials  by 
E.  Orowan  [^26^  and  G.  Irwin  [_2lJ'     This  work,  however, 
deals  entirely  with  nonimpulsive  loads,   and  is  unsuit- 
able for  application  to  the  present  case.     N.  F.  Mott 

made  the  initial  attempt  at  presenting  a  rational 
theory  for  the  fracture  of  cylindrical  ring  bombs.  He 
based  his  development  of  the  theory  upon  a  property  of 
the  materials  in  static  tension  tests,  known  as  "scatter". 
Apart  from  the  fact  that  materials  can  sustain  consider- 
ably larger  stresses  for  short  periods  of  time  than  under 
static  conditions  and  as  such  any  extension  of  static 
results  to  impulsive  situation  is  uncertain,  Mott's 
theory  is  dependent  upon  the  following  two  assumptions 
i.e.,    (1)   fracture  propagation  is  instantaneous  as 
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compared  to  rate  of  strain.    ^  Rate  of  strain  =  V/R, 
where  V  denotes  the  particle  velocity  and  R  is  the 
radius  of  the  cylinder  )    (2)  fracture  occurs  due  to 
exhaustion  of  ductility  (in  other  words,   fracture  is 
determined  by  total  strain  and  not  strain  rate) . 

An  alternative  theory  for  fracturing  of  cylindri- 
cal shells  is  presented  by  J.  S.  Rinehart  and  J.  Pearson 

Although  this  theory  agrees  with  that  of  Mott 
insofar  as  the  number  of  fragments  is  supposed  to  be 
proportional  to  radial  particle  velocity,   its  development 
is  based  upon  an  entirely  different  and  empirical  concept 
of  "critical  velocity  of  straining."     Velocity  of  strain-' 
ing  is  defined  as  the  rate  of  stretching  of  the  surface 
of  the  shell  and  is  equal  to  (27Tv).     The  critical  velocity 
of  straining  corresponds  to  the  minimum  value  of  V  for 
which  fracture  will  occur.     In  the  following,  their  basic 
formulation  shall  be  employed. 

At  the  outset,   it  is  useful  to  point  out  the 
distinctive  features  of  impulsive  loads.     Under  static 
loads,   stresses  and  strains  are  distributed  throughout 
the  body  and  local  stress-concentrations  are  important. 
However,  under  impulsive  loading  transient  stresses  and 
strains  may  exist  in  highly  localized  regions,   and  the 
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fracture  may  not  even  propagate  before  another  stress 
situation  comes  into  play. 

Impulsive  loads  set  up  stress  waves  in  the 
medium.     In  the  case  of  loading  such  as  that  occurring 
in  the  present  case  (blast  wave  hitting  head-on  into  the 
shell),   the  waves  generated  are  longitudinal.     If  the 
wave  is  one  of  compression,  the  particle  velocity  is  in 
the  same  direction  as  the  wave.     And,   if  the  wave  is  one 
of  tension,  the  particle  velocity  is  in  the  direction 
opposite  to  that  of  the  wave.     No  information  on  formulae 
for  wave  propagation  in  thin  spherical  shells  is  available. 
The  stress-wave  propagation  velocity  may  be  anywhere  be- 
tween 3,000  -  21.000  feet/seconds  (note  that  stress-wave 
propagation  velocity  is  not  the  same  as  particle  velocity) . 
This  should  be  compared  to  the  fact  that  velocity  of 
fracture  is  rarely  more  than  5,000  feet/seconds.     M.  Kolsky 

reported  experiments  on  polystyrene  plastic  rods 
in  which  fracture  starting  along  the  axis  did  not  extend 
to  the  surface.     Apart  from  the  fact  •of  stress-concentra- 
tion along  the  axis,   this  phenomenon  may  be  attributed 
to  the  rapidly  changing  stress-situation,   and  hence, 
strain- rate,     m  any  case,  it  is  clear  from  the  foregoing 
that  it  is  not  safe  to  assume  fracture  velocity  as  much 
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greater  than  velocity  of  straining  under  impulsive  load- 
ing.    This  is  especially  true  in  the  case  of  ductile 
materials  where  fracture  velocity  is  extremely  low. 

3 .     On  Shell  Fragmentation 

In  order  to  study  fracture  of  the  shell  it  is 
first  necessary  to  specify  the  loading.     The  shell  is 
supposed  to  be  loaded  by  detonating  a  chemical  explosive 
at  its  geometric  center.     Ignition  of  the  explosive  at 
the  geometric  center  of  the  shell  will  give  rise  to  a 
spherical  blast  wave.     Assuming  that  the  shell  has  no 
effect  on  the  detonation  wave,   the  shock-wave  collision 
with  the  shell  will  result  in  the  latter 's  being  loaded 
impulsively.     The  pressure  rises  instantaneously  to  its 
maximum  value  and  then  decays  to  zero  in  a  very  short  time. 
It  is  thus  clear  that  the  loading  here  is  not  only  dynamic 
but  also  impulsive. 

The  collision  of  the  shock  wave  with  the  shell  sets 
up  a  compression  wave  in  the  shell.     This  is  illustrated 
in  Figure  10a.     When  this  wave  reaches  the  outer  surface 
it  is  reflected  back  as  a  tension  wave  and  one  has  the 
situation  shown  in  Figure  10b.     The  particle  velocity  is 
now  obtained  by  superposition  and  is  2V .     Thus,  after 
the  shock  wave  strikes  the  shell  surface,  the  shell  is 
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Figure  10:     Propagation  of  Stress-Waves  in  Media 
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expanding  at  a  very  rapid  rate.     Following  Rinehart  and 
Pearson,  a  critical  velocity  of  straining  (V^^)  is 
postulated.     Fracture  will  occur  if  the  velocity  of 
straining  reaches  the  critical  velocity  of  straining.  If 
the  spherical  shell  is  radially  expanding  with  a  velocity 
V,  then  the  rate  of  increase  in  surface  area  is  given  by 
8  TT  R  V.     The  velocity  of  straining  is  now  defined  as 
equal  to  8  TT  V  and  the  critical  velocity  of  straining  as 
the  one  corresponding  to  the  minimum  value  of  V  for  which 
fracture  will  occur. 

The  question  that  must  be  answered  now  is  whether 
the  fracture  will  be  of  a  shear  (ductile)  or  tension 
(brittle)   type.     A  tension  fracture  usually  occurs  during 
expansion  of  the  shell.     A  shear  fracture,  on  the  other 
hand,  will  occur  during  release  of  compressive  stress  by 
the  reflected  tension  stress  wave.     If  fracture  occurs 
during  first  compression  wave,  it  will  be  of  a  brittle 
type.     If,  however,   fracture  does  not  occur  during  the 
first  compression  wave,   it  may  turn  out  to  be  a  compli- 
cated mixture  of  brittle  and  ductile  types.     This  latter 


type  has  been  observed  in  thick  cylindrical  shells  29 


The  fracture  behavior  may  also  be  further  complicated  by 
relatively  high  temperatures  reached  at  the  shock  front 
[^6^  and  consequent  heating  of  the  shell  surface. 
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4 .     Experimental  Program 

In  light  of  Rinehart  and  Pearson's  expression  for 
cylindrical  shells,   the  following  relationship  is  pro- 
posed for  the  number  of  fragments: 


A  number  of  experiments  were  conducted  on  polystyrene 
shells  of  several  diameters.     The  shell- thickness  was 
approximately  0.05  inches.     Shells  were  loaded  internally 
by  detonating  a  PBRDX  explosive.     The  combination  used 
prior  to  assembly  is  reproduced  in  Figure  11.     Two  differ- 
ent sizes  of  charges  were  employed  in  the  testing  program 
i.e.,    (1)  a  spherical  charge  of  3/8  inch  diameter  varying 
in  weight  from  0.56  grams  to  0.62  grams  and  (2)  a  spheri- 
cal charge  of  1/2  inch  diameter  varying  in  weight  from 
1.27  grams  to  1.31  grams.     These  charges  were  mounted  at 
the  geometric  center  of  the  spheres  when  assembled  from 
two  matching  hemispheres. 

The  experimental  program  was  conducted  in  two  parts. 
The  first  part  consisted  of  experimental  firings  for  high- 


(40) 


A  logical  experimental  program  is  now  to  determine 


the  ratio 


Figure  11:     Combination  of  Shell,  Detonator 
and  Charge  Prior  to  Assembly 
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speed  measurements  of  case  break-up.     The  second  part 
consisted  of  similar  experimental  firings  for  recovery  of 
fragments.     It  was  not  possible  to  obtain  pictures  and  to 
recover  fragments  from  the  same  experiment  because  of  the 
need  for  an  intense  light  source.     The  two  sections  are 
described  below. 

The  photos  were  taken  with  a  Beckman  and  Whitley 
framing  camera,  model  189  operating  at  500,000  frames 
per  second.     The  duration  of  the  individual  exposure  was 
approximately  0.3  microseconds.    Photo  illumination  was 
provided  from  a  f ourteen-inch  cone  filled  with  argon  gas 
and  initiated  by  a  100  gram  charge  of  composition  C-3. 
A  typical  photo- sequence  is  reproduced  in  Figure  12. 

The  recovery  of  fragments  was  never  100  per  cent  as 
the  experimental  firing  bay  is  open  on  one  side.  However, 
it  would  appear  that  only  a  very  small  amount  (less  than 
3-5  per  cent)  was  lost  as  the  weight  of  recovered  frag- 
ments in  the  case  when  only  the  detonator  was  used 
equalled  the  weight  prior  to  assembly.     The  discrepency 
between  the  weight  before  firing  and  the  weight  of  the 
recovered  fragments  is  because  of  the  evaporation  and 
burning  of  the  plastic  as  a  result  of  the  intense  heats 
reached  at  the  shock-front  hitting  the  shell.     It  is  to 
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4t 

Figure  12:     Fracture  of  1  11/16  Inch  Diameter  Plastic 
Sphere.     Charge  Weight:     1.3  grams. 
Frame  Rate:     500,000  Frames/Second 
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be  observed  that  nearly  one  third  of  the  plastic  weight 
either  evaporated  or  burned  in  each  case.     The  fragments 
recovered  were  generally  of  three  sizes  i.e.,    (i)  very 
large,    (ii)  moderate  and  (iii)  extremely  small  bordering 
on  fine  powder.     The  number  of  fragments  listed  in 
Table  6  is  nominal.     It  was  found  convenient  to  count  the 
number  of  largest  fragments  equalling  two  thirds  of  the 
recovered  weight  and  then  multiply  the  number  by  3/2. 
In  justification  of  this  process,  it  may  be  said  that 
very  small  fragments  may  be  neglected  in  first  approxima- 
tions.    The  value  of  particle  velocity  V  was  obtained 
directly  from  high-speed  photographs.     The  last  two  columns 
in  Table  6  represent  calculated  values. 

The  purpose  of  the  experiments  was  two- fold  i.e., 
(1)  to  study  the  nature  of  fracture  and  (2)  to  determine 
A*.     The  plastic  spheres  exhibited  a  great  isotropy  in  ' 
fracture  (see  Figure  12)  revealing  a  predominantly  brittle 
fracture.     The  value  of  ^  *  was  calculated  from  observed 
values  of  V  and  N.     It  is  seen  from  Table  6  that  this 
value  is  the  same  within  experimental  limits,   for  the 
shell  of  a  particular  diameter  regardless  of  the  charge 
size.     This  is  as  it  should  be. 

A  result  of  still  greater  importance  is  the  fact  that 
D    ^*  is  a  constant.     Table  6  reveals  that         ^  *  varies 
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by  less  than  10  per  cent  around  a  median  value  of  0.549. 
This  agreement  is  extremely  good  in  view  of  the  various 
experimental  limitations,  such  as,    (1)  difference  in 
charge-sizes,    (2)  failure  to  center  charges  properly, 
(3)  necessity  of  obtaining  pictures  and  fragments  from 
separate  experiments,    (4)  failure  to  collect  all  frag- 
ments and  (5)  error  in  reading  high-speed  photographs. 
It  goes  without  saying  that  elimination  of  one  or  several 
of  these  sources  of  error  should  yield  substantially  better 

agreement  in  the  value  of  ,     it  follows  then  from 

3 

constancy  of  D    /^*  that  the  critical  velocity  of  strain- 
ing (V^j.)  is  directly  proportional  to  the  cube  of  the 
shell  diameter. 

A  part  of  the  experimental  program  reported  in 
Table  6  was  repeated  by  filling  the  empty  space  between 
the  charge  and  the  plastic  shell  with  water.     The  results 
are  reproduced  in  Table  7.     They  show  a  remarkable  agree- 
ment With  the  earlier  results,  confirming  thereby  once 
again  the  hypothesis  that  /\^*  is  independent  of  pressure 
for  a  particular  shell-size.     The  present  tests  show  a 
slightly  higher  value  for        X*  •     The  difference  is  not, 
however,  significant  as  the  upper  limit  in  the  former 
case  and  the  lower  limit  in  the  present  case  cross  each 
other. 
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5 .     Some  Dimensional  Considerations 

In  the  preceding  section,  It  was  shovm  that  V 

cr 

is  directly  proportional  to  D^.     Now  the  effect  of  wall 

thickness  (W. )and  material  properties  on  V      shall  be 
t  •*■  cr 

examined  on  the  basis  of  dimensional  analysis. 

Since  the  shells  under  consideration  are  thin,  shear 
effects  may  be  neglected.     The  primary  dimensional  parame- 
ters are  then  (1)   shell  diameter  (D) ,    (2)  wall  thickness 
(W^) /    (3)  Young's  modulus   (E) ,   and  (4)  density  of  material 
(  ^) .     The  dependence  of  V^^  on  these  parameters  may  be 
expressed  mathematically  as 

(41) 

Assuming  that  it  is  possible  to  expand  f  in  powers  of 
W^,  E  and    ^ ^,  (41)  may  be  written  as: 


where   X^,  >  ,  2^'   '  V     are  nondimensional  constants  and 


(42) 


Then  from  dimensional  analysis 
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Where  M,  L,  T  respectively  denote  dimensions  of  mass, 
length  and  time.     Comparing  exponents  of  M,  L  and  T  on 
both  sides,  one  finds 

V'=  -3 
-O"  V2 

Therefore  from  equations   (41)  and  (42)  there  follows 


m 

(43) 

It  is  now  interesting  to  examine  how  the  variables 
(D,  W^,  E,  effect  the  number  of  fragments.     For  a  speci- 

fied impulse,  particle  velocity  V  is  inversely  proportional 
to        (M  is  mass  of  the  shell)  and  hence  to  W^^.  However, 

it  is  independent  of  material  properties  E  and   5'  .  Then 

m* 

from  (40) ,   it  becomes  obvious  that  wall  thickness  has  no 
effect  on  the  number  of  fragments.     Also,   if  one  denotes  by 
siabscripts  1  and  2  two  different  materials,  then 


(44) 

Although  the  above  results  are  derived  from  quite  general 
considerations,  it  must  be  emphasized  that  they  should  be 
used  with  caution  pending  experimental  confirmation  as 
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several  nondimensional  parameters  (such  as  Polsson's  ratio) 
may  also  effect  the  behavior  of  shells, 

6.  Conclusions 

Experimental  tests  on  polystyrene  shells  confirm 
the  possibility  of  extending  the  eir^irical  criteria  of 
Rinehart  and  Pearson  for  cylindrical  shells  to  the  case  of 
spherical  shells.     Experimental  results  also  prove  that 
^  *  is  dependent  only  upon  shell-size  and  material  and  not 
upon  the  charge-size.     Further,  it  was  observed  that  V 

cr 

is  directly  proportional  to  the  cube  of  the  shell-diameter. 
The  fracture  pattern  was  found  to  be  isotropic  and  brittle. 
It  is  shown  from  dimensional  analysis  that  wall  thickness 
has  no  effect  on  the  number  of  fragments.     The  effect  of 
material  properties  is  also  investigated. 


CHAPTER  IV 

ON  THE  CALCULATION  OF  FRACTURE-TIME  FOR  THIN 
SPHERICAL  SHELLS  SUBJECT  TO 
INTERNAL  BLAST  LOADING 

1.  Introduction 

In  Chapter  I,  a  similarity  solution  for  obtaining 
the  pressure  behind  the  shock  front  in  a  spherically 
synunetric  blast  wave  generated  by  a  finite  charge  was 
outlined.     This  method  was  used  in  conjunction  with 
experimental  data  to  find  the  pressure  distribution  in  the 
explosive  gases  behind  the  contact  front  in  Chapter  II". 
In  Chapter  III,  the  nature  of  fracture  under  impulsive 
loading  was  examined.     An  empirical  formula  for  the  number 
of  fragments  was  proposed.     It  was  shown  that  the  number 
of  fragments  (N)  was  directly  proportional  to  the  particle 
velocity  V  (P)  where  P  is  the  dynamic  pressure,     in  the 
present  Chapter,   the  relationship  between  V  and  P  is 
examined  through  the  introduction  of  the  concept  of 
fracture  time.     Fracture  time  is  defined  as  the  time  which 
must  elapse  between  the  impact  of  the  blast  wave  on  the 
shell  and  the  rupture  of  the  shell.     The  general  relation 
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between  V  and  P  is  simplified  by  introduction  of  thin-shell 
assumptions.     Finally,   the  relation  is  applied  to  the 
theoretical  and  experimental  data  previously  reported  in 
Chapters  I  -  III. 

2.     Preliminary  Considerations 

The  blast  wave  on  hitting  the  shell  imparts  to  it  an 
impulse.     This  can  be  easily  calculated  by  integrating 
[pit)  +   "5*  (t)  u2(t)J  from  the  time  of  impact   (t=0)  to  the 
time  of  rupture   (t  )       3lJ.  Mathematically; 

^  ^  (45) 

where  R  is  the  radius  of  the  shell,  p(t)  the  pressure 
acting  on  the  surface  of  the  shell,  ^  (t)   the  density  of 
the  impinging  gases  and  u(t)  the  velocity  of  the  gases. 

The  impulse  imported  to  the  shell  appears  as  the 
momentum  (MV)  of  the  shell.  Thus 


(46) 

with 


M  denotes  the  mass  of  the  shell  and  V  the  particle  velocity. 
Relation  (46)  is  quite  general.     No  assumptions  have  been 
introduced  except  for  the  fact  that  in  proposing  equation 
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(46)  it  is  implicitly  recognized  that  the  energy  dissi- 
pated in  fracture  is  negligible  in  comparison  to  the 
total  energy  associated  with  the  impulse.     At  this  stage, 
certain  assumptions  are  introduced  in  order  to  simplify 
equation  (46) .     The  first  two  are:     (i)  The  blast  wave 
is  not  affected  by  the  shell  i.e.  reflection  may  be 
ignored.     (ii)  The  dynamic  pressure  P  remains  constant 
during  the  fracture  time  t^. 

It  is  hard  to  justify  these  assumptions  purely  from 
theoretical  reasoning.     Later  on,  more  will  be  said  on 
this  point.    For  the  present,  this  shall  be  accepted  as 
a  working  hypothesis. 

So  far,  nothing  was  said  about  the  shell  thickness. 
In  the  following,  under  thin-shell  assumptions,  the 
structure  of  equation  (46)  will  be  examined, 

3.     Thin-Shell  Simplifications 

The  impact  of  the  blast  wave  on  the  shell-surface 
sets  up  longitudinal  stress  waves  in  the  shell.  The 
stress-wave  propagation  velocity  may  lie  anywhere  between 
3,000  -  21,000  feet/second.       Thus,   the  time  for  the  head 
of  the  stress-wave  to  travel  a  distance  of  0.05  inch 
(typical  shell  thickness  employed  in  experiments  reported 
in  Chapter  III)  is  of  the  order  of  (0.2  -  1.4) ^sec.  It 
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was  observed  that  the  fracture  of  the  shell  was  primarily 
brittle  and  that  it  exhibited  no  scabbing.     This  implies 
that  the  fracture  occurred  during  the  onward  first 
compression  wave  and  before  its  reflection  as  a  tension 
wave  at  the  outer  surface  of  the  shell  [29].     Hence,  the 
fracture  time  (tf)  must  lie  in  the  range  (0.2  -  1.4) 
^sec.     For  this  order  of  time,  no  great  decay  occurs  at 
the  shock  front  (Chapter  II) .     The  decay  is  also  balanced 
by  reflection  of  the  shock  wave  from  the  surface  of  the 
shell.     It  would  seem  reasonable,   therefore,   to  take 
values  of  p  as  that  existing  just  behind  the  shock  front. 
The  values  of    ^  and  u  are  taken  as  those  existing  across 
the  contact  front  as  gases  in  the  contact  front  have  a 
much  higher  density  than  compressed  air.     Then  equation 
(46)  becomes 

and,  hence, 

Ml/ 


(47) 


4.     Experimental  Results 

Experimental  results,   involving  the  measurement  and 
calculation  of  the  pressure  at  the  shock  front  resulting 
from  exploding  PBRDX  charges,  were  reported  in  Chapter  II, 
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The  first  six  columns  in  Table  8  are  based  upon  these 

* 

results.      Columns  7  and  8  are  reproduced  from  Chapter  III 
(Table  6) ,  where  the  behavior  of  a  plastic  shell  under 
blast  loading  (charges  employed  were  identical  with 
those  in  Chapter  II)   is  discussed.     The  last  column 
reports  fracture  time  calculated  by  employing  equation 
(47). 

It  is  seen  that  the  fracture  time  t^  is  consistently 
less  for  the  larger  charge  indicating  that  the  fracture 
time  is  dependent  upon  dynamic  pressure  P.     This  is  to 
be  expected  as  the  impulse  I  imparted  to  the  shell  and 
the  particle  velocity  V  depend  upon  the  dynamic  pressure 
P. 

The  fracture  time  is  of  the  order  of  (0.6-2.0)  ^isec. 
This  is  approximately  twice  what  one  would  expect  from 
theoretical  analysis  of  stress  wave  propagation.  The 
discrepancy  is  due  to  the  neglect  of  reflection.     If  even 
a  simple  reflection  law  is  accepted,   it  would  double  the 
dynamic  pressure   (P)  and  thus  reduce  the  fracture  time 
by  half. 


In  Table  8,   all  the  units  employed  belong  to  the 
CGS  system  regardless  of  units  used  in  Chapters  II  and  III. 
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CHAPTER  V 
CONCLUSIONS 

On  the  basis  of  preceding  mathematical  and  empiri- 
cal analysis,  following  conclusions  may  be  drawn: 

(1)  There  exists  a  group  of  similarity  transforma- 
tions for  the  case  of  axisymmetric  flow  when  the  Euler 
equations  are  given  in  spherical  coordinates.     It  is 
possible  to  employ  this  set  of  transformations  along 
with  empirical  data  to  obtain  approximate  solutions  for 
the  spherically  symmetric  and  axisymmetric  blast  waves 
generated  by  point-source  or  finite-source  explosions. 

(2)  Dominant  fracture  pattern  in  thin  polystyrene 

spherical  shells  under  internal  blast  loading  is  brittle. 

The  number  of  fragments   (N)  is  inversely  proportional  to  t 

critical  velocity  of  straining  (V    ).     It  is,  however, 

cr 

independent  of  the  wall  thickness  of  the  shell.  Also, 
critical  velocity  of  straining  (V^j.)  is  directly  propor- 
tional to  the  shell-diameter  (D) . 

(3)  The  fracture- time  (t^)  depends  upon  dynamic 
pressure  (P)  and  decreases  with  increasing  P. 
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